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ABSTRACT 

Observational and theoretical evidence suggests that high-energy Galactic cosmic rays are primarily 
accelerated by supernova remnants. If also true for low-energy cosmic rays, the ionization rate near 
a supernova remnant should be higher than in the general Galactic interstellar medium (ISM). We 
have searched for H|^ absorption features in 6 sight lines which pass through molecular material 
near IC 443 — a well-studied case of a supernova remnant interacting with its surrounding molecular 
material — for the purpose of inferring the cosmic-ray ionization rate in the region. In 2 of the sight 
lines (toward ALS 8828 and HD 254577) we find large H;^ column densities, 7V(H;^) w 3 x lO" cm'^, 
and deduce ionization rates of ^2 ~ 2 x lO"""^^ s~^, about 5 times larger than inferred toward average 
diffuse molecular cloud sight lines. However, the 3cr upper limits found for the other 4 sight lines are 
consistent with typical Galactic values. This wide range of ionization rates is likely the result of particle 
acceleration and propagation effects, which predict that the cosmic-ray spectrum and thus ionization 
rate should vary in and around the remnant. While we cannot determine if the HJ absorption arises 
in post-shock (interior) or pre-shock (exterior) gas, the large inferred ionization rates suggest that IC 
443 is in fact accelerating a large population of low-energy cosmic rays. Still, it is unclear whether 
this population can propagate far enough into the ISM to account for the ionization rate inferred in 
diffuse Galactic sight lines. 

Subject headings: astrochemistry - cosmic rays - ISM: supernova remnants 



1. INTRODUCTION 

As cosmic rays propagate through the interstellar 
medium (ISM) they interact with the ambient material. 
These interactions include excitation and ionization of 
atoms and molecules, spallation of nuclei, excitation of 
nuclear states, and the production of neutral pions (tt") 
which decay into gamma-rays. Evidence suggests that 
Galactic cosmic rays are primarily accelerated by super- 
nova remnants (SNRs) through the process of diffusive 
shock acceleration fe.g. lDrurvlll983tli31andford fc Eichleri 
Il987f) , so interstellar clouds in close proximity to an SNR 
should provide a prime "laboratory" for studying these 
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interactions. IC 443 represents such a case, as portions 
of the SNR shock are known to be interacting with the 
neighboring molecular clouds. 

I C 443 is an inte rmediate age remnant (about 30,000 
yr; iChevalieTI Il999t ) located in the Galactic anti-center 
region (I, h) « (189°, -1-3°) at a di stance of about 1.5 kpc 
in the Gem OBI association (Wel sh fc SallmenI [2003). 
and is a particularly well-studied SNR. Figure [1] shows 
the red image of IC 443 taken during the Second Palo- 
mar Observatory Sky Survey. The remnant is composed 
of subshells A and B; shell A is to the NE — its cen- 
ter at a = 06^17'"08.4", 5 = -|-22°36'39.4" J2000.0 is 
marked by the cross — while shell B is to the SW. Adopt- 
ing a distance of 1.5 kpc, the radii of subshells A and 
B are about 7 pc and 11 pc, respectively. Between the 
subshells is a darker lane that runs across the remnant 
from the NW to SE. This is a mo lecular cloud which 
has been mapped in ^^CO emissio n (iCor nett et aLllTOTTl : 
iDickman et al.|[T992l : iZhang et al. 2009) , and is known to 
be in the foreground because it absorbs X-rays emitted 
by the hot remnant interior (|Troia et al.l 120061 ). Aside 
from this quiescent foreground cloud, observations of the 
J = 1 — line of ^^CO also show s hocked molec- 
ular material coincident with IC 443 (DeNovcr 19791 
iHuang et aITl986;.Dickman et al.ll992; Wang k Scovillg 
1992) . These shoc ked molecular clumps fi rst identified 
bv iDeNoveil (|1979D and IHuang et al.l (|1986D in CO have 
also been observ ed in several atom i c and small molec- 
ular specie s (e.g. [W hite et al.' '1987'; Burton et al. 1988t 
Ean Dishoeck_et3!fll993; White 1994; Snell et al. 2005), 
and are thought to be the result of the expanding SNR 
interacting with the surrounding ISM. While many of the 
shocked clumps are coincident with the quiescent gas, it 
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is unclear whether or not they are part of the foreground 
cloud (i.e. the back portions of the foreground cloud are 
beginning to interact with the SNR blast wave) , or if the 
foreground cloud is separated from IC 443. 

Chemical analyses performed in various studies of the 
shocked clumps around IC 443 su ggest a n enha nced ion- 
ization rate due to cosmic rays. iWhitd (|1994D found a 
C/CO ratio much higher than in typical dense clouds 
and concluded that shocks and/or a la rge flux of cos- 
mic rays m ust be responsib le. Both iClaussen et al.l 
(jl99^ and H ewitt eraLl (pfiOfih observed OH (1720 MHz) 
masers toward some of these clumps. It is thought 
that this OH is formed when the free electrons pro- 
duced during ionization events collide with and excite 
H2, which in turn em its UV photons that dissociate H2O 
(jWardle fc Yusef-^Za dch 2002). In order to convert nearly 
all of the H2O into OH, thus generating the large column 
of OH necessary to produce the observed masers, a high 
ionization rate due to X-rays and/or cosmic rays is re- 
quired. Estimates of the ionization rate due to X-rays 
( Yusef-Zadeh et al.|[20Ql and cosmic rays ()Hewitt et al.l 
|2009[) near IC 443 are similar (a few times 10-^^ s'^), so 
it may be that both play a role in generating OH. How- 
ever, none of these analyses alone can determine exactly 
how important cosmic-ray ionization and excitation are 
to the processes considered. 

Recently, many studies of IC 443 have focused on 
the production of pionic gamma-rays via interactions 
between hadronic cosmic rays and ambient nucleons. 
Gamma-ray observations of IC 4 43 have be en performed 
by E GRET (lEsposito et allHoOfih . MA GI G (lAlbert et al.l 
mMi. VERITAS (lAcciari et al.l [2009ll. FermA LAT 
(|Abdo et alllMoh . and AGILE (iTavani et al. l [20Tflh . All 
show gamma-ray emission that appears to be coincident 
with gas in the vicinity of IC 443, thus supporting an en- 
hanced cosmic ray flux in the region. Because produc- 
tion requires cosmic-ray protons with Ekin > 280 MeV, 
gamma-ray observations cannot constrain the cosmic-ray 
flux at lower energies. 

To investigate the flux of lower-energy cosmic rays, 
we study the cosmic-ray ionization of H2, a process 
dominated by protons with 1 MeV < Ekin < 1 GeV 
(jlndriolo et al.ll200l IPadovani et al.l[2009( ). The ioniza- 
tion rate of H2, (^2, can be inferred from observations 
of H;^ assuming a rather simple chemical network. H2 is 
first ionized, after which the ion collides with another H2, 
thus forming H;^. Either dissociative recombination with 
electrons (diffuse clouds) or proton transfer to CO, O, 
and C (dense clouds) are the primary destruction routes 
for H;^ depending on the environment. In this paper, we 
present observations searching for absorption lines of 
along sight lines which pass through molecular material 
near IC 443. We then use the results of these observa- 
tions in combination with the simple chemical scheme 
outlined above to infer the cosmic-ray ionization rate of 
H2. 

2. OBSERVATIONS 

This project examined 6 target sight lines toward the 
stars ALS 8828, HD 254577, HD 254755, HD 43582, HD 
43703, and HD 43907, all of which are shown in Figure 
[1] to the immediate left of the labels A-F, respectively. 
Target selection was based on various criteria, includ- 



ing L-band magnitude, previously detected molecules, 
and evidence t hat the background stars were in fact be- 
hind the SNR (|Welsh fc SallmenI [20031: IHirschauer eFall 
120091 ). Basic properties of the se sight lines are avail- 
able in IHirschauer et al.l (|2009( ). Observations focused 
primarily on transitions arising from the (J, K) = (1, 1) 
and (1,0) levels of the ground vibrational state of H;^, 
the only levels significantly populated at average diffuse 
cloud temperatures (T ~ 60 K). Transitions from higher 
energy levels (e.g. (2, 1) and (3, 3)) were covered as al- 
lowed by the instrument, but absorption at these wave- 
lengths was not expected. 

Spectra were obtained us hig the Near Infrared Echelle 
Spectrograph (NIRSPEC; iMcLean et all I1998D at the 
W. M. Keck Observatory, and the Infrared Camera and 
Spectrograph (IRCS; Kobayashi et al. 2000) at the Sub- 
aru Telescope. All NIRSPEC observations were per- 
formed on 2009 Nov 5 and 6 using the 3-pixel (0.432" ) slit 
to provide a resolving power of ^ 25, 000. The KL filter 
was used in combination with the echelle grating/cross- 
disperser settings of 64.82/33.5 in order to simultane- 
ously cover the R{1, 1)", R{1, 0), and -R(3, 3)' transitions. 
All IRCS observations were performed on 2009 Dec 12 
and 13 in echelle mode using the 2-pixel (0.14") slit to 
provide a resolving power of ~ 17, 300. The adaptive op- 
tics system (A0188) was utilized in order to maximize 
starlight passing through the narrow slit. The L filter 
was used in combination with the echelle grating/cross- 
disperser settings of 8350/6100 in order to simultane- 
ously cover the i?(l,l)", i?(l,0), -R(l, 1)', (5(1,1), and 
(5(1,0) transitions. Select properties of the six targeted 
transitions are listed in Table [U and their locations with 
respect to atmospheric absorption features are shown in 
Figure [2] A log containing the list of science targets 
and exposure times for each night is shown in Table [2j 
In addition to the science targets, the bright, early-type 
stars (3 Tau and C Tau were observed for use as telluric 
standards. These standards were observed so that they 
would have good air mass matches with the science tar- 
gets, and to account for atmospheric variability over the 
course of each night. For all observations, the star was 
nodded along the slit in an ABBA pattern in order to 
facilitate the removal of atmospheric emission lines and 
dark current via the subtraction of neighboring images. 

3. DATA REDUCTION 

Our data reduction process combines the use of stan- 
dard IRAlO procedures and macros written in IGOR 
PrcO. Due to differences between the NIRSPEC and 
IRCS observations, each data set required slightly differ- 
ent reduction techniques. 

3.1. Keck Data Reduction Process 

A bad pixel map was created from the average of sev- 
eral dark frames, and these pixels were interpolated over 
in the flat-field and object images. Images were then cut 
into two sections, with each section containing one of the 
orders of interest from the cross-dispersed spectrograph. 
These sections were then treated as individual images for 
the remainder of the reduction process. The flats were 

http;//iraf.noao.edu/ 
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then combined, and each object frame was divided by 
the normalized, averaged flat-field. Neighboring AB im- 
age pairs were subtracted from each other to remove at- 
mospheric emission and dark current. One-dimensional 
spectra were extracted for each order using apall, and 
imported to IGOR Pro. 

3.2. Subaru Data Reduction Process 

While no dark frames were taken at Subaru, calibra- 
tion frames with the lamp off were used to create a bad 
pixel map. These pixels were interpolated over in the 
object and flat-field frames, and again each image was 
cut into two sections containing an order of interest. In- 
spection of the averaged, normalized flat-field frame for 
each order showed a low signal-to-noise ratio (S/N), so 
object frames were not divided by these flats. As before, 
one-dimensional spectra were extracted with apall and 
imported to IGOR Pro. In order to remove a saw-tooth 
pattern (the result of different readout channels) from 
these spectra, a moving average was taken for both the 
odd and even numbered pixels, and the even pixels were 
then scaled by the ratio of these averages. 

3.3. Shared Reduction Processes 

Individual spectra within an exposure sequence for a 
given target were then added together. In this pro- 
cess, spectra with S/N much lower than average (due 
to cirrus clouds or a bad nodding sequence) were ex- 
cluded. Each summed spectrum was then divided by a 
telluric standard to remove atmospheric absorption fea- 
tures and normalize the spectrum. These ratioed spec- 
tra were wavelength calibrated with a typical accuracy of 
~ 2 km using the vacuum wavelengths of the atmo- 
spheric absorption lines. Calibrated spectra were then 
shifted into the local standard of rest (LSR) frame. At 
this point all spectra of a given target from a single tele- 
scope were combined using a variance-weighted mearP^ 
= Er=o(^i/^i^)/ Er=o(l/^i^)' where ctj is the stan- 
dard deviation on the continuum near the H;^ line posi- 
tions) to produce a final Keck and Subaru spectrum for 
each sight line. The i?(l,l)', Q(l,l), and Q(1,0) tran- 
sitions were only covered at Subaru, while the i?(3,3)' 
transition was only covered at Keck, so these are the final 
spectra presented in FigurelH The i?(l, I)" and i?(l,0) 
transitions were covered at both telescopes, so the final 
Keck and Subaru spectra can be combined to obtain a 
higher S/N. The Keck spectra were interpolated onto the 
lower resolution Subaru wavelength scale, and all spectra 
were again combined via the weighting scheme described 
above. The resulting spectra for our six target sight lines 
are shown in Figure |31 

4. RESULTS 

It is clear from Figure |3] that absorption is only 
detected in the sight lines toward ALS 8828 and HD 
254577. The i?(l,l)" and i?(l,0) lines are quite strong 
toward ALS 8828. Absorption from the higher energy 
(3,3) and (2,1) states was not detected, as expected given 
diffuse molecular cloud conditions. The sight line toward 
HD 254577 shows absorption from the i?(l, 1)", i?(l, 0), 

This averaging scheme is equivalent to weighting each spec- 
trum by (S/N)2. 



i?(l,l)' and Q(1,0) transitions of H;^. Absorption due 
to the (5(1, 1) transition must also be present, but it is 
not detected. This is probably the result of 3 factors: 
(1) the intrinsic strength of the (5(1, 1) transition is the 
weakest of the 5 transitions examined (see Table [T]); (2) 
the spectrum near the (5(1, 1) transition has a low S/N 
due to lower illumination of the echelle order in which 
it appears; (3) the (5(1, 1) transition is overlapped by a 
strong atmospheric N2O line (see panel 4 of Figure [2]), 
making removal of telluric features uncertain. Imperfect 
removal of this atmospheric line is also the most likely 
cause of the feature in the (5(1, 1) spectrum of HD 254755 
that appears at the expected velocity. This feature can- 
not be due to Hg , as there is no absorption by any of the 
other stronger H|^ transitions which arise from the same 
state. The positive spike near 130 km s""'^ in the i?(l, 1)' 
spectra of HD 254577 and HD 254755 is an instrumental 
artifact. Spectra of HD 43582, HD 43703, and HD 43907 
also show no absorption features from H^. 

5. ANALYSIS 

Equivalent widths were determined using Gaussian fits 
to the absorption features. Uncertainties were deter- 
mined from the standard deviation, a, on the residual 
continuum after subtracting by the Gaussian line pro- 
files. Interstellar gas velocities and velocity full width at 
half-maxima (FWHM) were also determined during this 
fitting procedure. In the case of non-detections, upper 
limits were determined from 3cr on the continuum across 
the expected position of a line assuming a FWHM of 
16 km s~^ (the resolution obtained with IRCS at Sub- 
aru). Column densities were derived from equivalent 
widths using the standard relation given optically thin 
absorption lines and the transition dipole moments and 
wavelengths listed in Table [TJ All of these results are 
shown in Table [3] 

Following the analysis of llndriolo et al.l (|2007| ). we 
adopt the simple chemical scheme in diffuse clouds where 
every ionization of an H2 molecule leads to H;^ , and dis- 
sociative recombination with electrons is the dominant 
mechanism by which H^ is destroyed. This results in 
the steady-state equation 

C2n(H2) = k,n,n{R+), (1) 

(|Geballe et al.l [19991 ) . where (2 is the ionization rate of 
H2, and fee is the electron recombination rate coefficient 
of H|^ . Substituting the electron fraction (defined as Xe = 
ne/nn, where riu = n(H)-f 2n(H2)) into equation ([T]) and 
solving for the ionization rate gives 

"(Hj) 

C2 = fceXenn-TTJ^. 2 

"-(II2) 

Although it would be desirable to trace the ionization 
rate as a function of position throughout the cloud, vari- 
ations in density along the line of sight cannot be de- 
termined via observations. Instead, we infer the aver- 
age ionization rate in a cloud by using average number 
densities. By definition, {n{ll^)) and (n(H2)) can be re- 
placed with N(R'^)/L and A^(H2)/L, respectively (where 
L is the cloud path length) , thus putting equation ^ in 
terms of observables. As Hj^ will form wherever there is 
an appreciable amount of II2 , it is reasonable to assume 
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that the path length for both species is the same, such 
that 



C2 = keXeU-n 



7V(H2) 



(3) 



Because the ratio 7i(H^)/n(H2) is not expected to 
vary widely in mod els of diffuse molecular clouds (e.g. 
INeufeld et aLl[2n05h . this should give a representative 
value of the ionization rate throughout the entire cloud. 

Assuming that the vast majority of electrons in dif- 
fuse molecular clouds come from photoionized carbon, 
can be approximated by N{C^)/N^ {N-^ is the col- 
umn density analog to nn), which was found to be 
about 1.5 X 10""^ .^^Q^^. multiple di f jFuse cloud sight 
lines (ICardelli et al.l IToM ISofia et al.ll200l . The hy- 
drogen number density can be estimated by both a 
rotation-excitation analysis of C2 observations and a 
restricted chemical ana lysis based on CN observations 
(|Hirschauer et alll2009D . The C2 analysis also gives a 
best-fit kinetic temperature which we use in calculat- 
ing fee. The temperature dependency of ke as deter- 
mincd from laboratory work is reported in IMcCall et al.l 
p004) . While molecular hydrogen has not been observed 
in absorption along any of our target sight lines, abun- 
dances of H 2 and CH tend to be linearly related in dif- 
fuse clouds HFedermanlllQSl [Mattila 1986; Sheffcr__eLal] 
[2OO8I) . We use the relationship derived from the largest, 
most recent data se t— jV(CH)/iV(H2) = 3.5^1,1 x 10"^ 
(jSheffer et al.l I2008[ ) — in combination with CH column 
densities reported bv lHirschauer et aLl (j2009D to estimate 
A^(H2). Finally, the total Tl^ column density is deter- 
mined by adding iV(l, 0) and iV(l, 1). These input values 
and/or the parameters on which they depend, as well as 
the inferred ionization rates are shown in Table 21 

While H2 is ionized by both cosmic rays and 
X-rays, most of the X-ray flux should be attenu- 
ated in a relatively thin layer at the cloud exterior 
(jGlassgold fc Langen 119741 ). The ionization rate due 
to X-rays at the ed ge of IC 443 was estini ated to be 
Cx = 3.6 X 10-16 s-i (|Yusef-Zadeh et al.|[200l . and must 
be much lower in cloud interiors. As a result, the ion- 
ization rates we infer should be primarily due to cosmic 
rays. 

6. DISCUSSION 

Having computed the cosmic-ray ionization rate for 
clouds in the vicinity of IC 443, we compare our results 
to those from previous studies. The aver age ionization 
rate in diffuse molecular clouds found by llndriolo et al.l 
(|2007D using H3 was C2 = 4 x IQ-^^ s-^, several times 
lower than found toward ALS 8828 and HD 254577 



(C2 = 16+^2 X 10"^^ and 26^11 x 10"^*^ s'^, respec- 
tively). In fact, the ionization rates inferred for these 2 
sight lines are more than twice the highest rates previ- 
ously found in diffuse molecular clouds toward ( Per and 
X Per (C2 « 7 X 10'^^ s'^). While ALS 8828 and HD 
254577 present exceptionally high ionization rates, the 
other 3 sight lines observed near IC 443 do not (due to 
the low S/N obtained toward HD 43907, the derived up- 
per limit for that particular sight line is not exceptionally 
meaningful, and so we exclude it from further consider- 
ation). Instead, the 3(t upper limits for C2 presented in 
Table|4]are consistent with ionization rates of a few times 



10^16 s^-^, typical of diffuse molecular clouds. These dif- 
ferences are quite striking, and warrant discussion. 

There are two plausible explanations for why Hj]" would 
be detected toward ALS 8828 and HD 254577 but not HD 
254755, HD 43582, HD 43703, and HD 43907, and they 
can most easily be seen when equation ([3]) is rearranged 
to show that A^(Hj^) oc C2/{xen-H)- Given this scaling we 
can posit that either the product Xgnn (i-e. the electron 
density) is lower along these 2 sight lines, or (2 is higher, 
and we examine these possibilities in turn. 

6.1. Lower Electron Density 

As stated in Section 5, we have assumed an electron 
fraction that is consistent with observations of C"*" in 
several diffuse molecular clouds. In denser environments 
though, the predominant form of carbon shifts from C^ 
to C, and eventually to CO, thus decreasing the electron 
density. Adopting a reduced value for Xenu requires a 
corresponding decrease in ^2 to match the observed H^ 
column density. It could then be argued that the en- 
hanced ionization rate we calculate for the 2 sight lines 
where we detect H|^ is actually just an artifact of not 
recognizing a decreased destruction rate. 

However, there are observations which seem to ar- 
gue against this possibility. The C2 rotation-excitation 
and CN restric ted chemical analyses performed by 
IHirschauer et al.l ([2)09) suggest densities of 200 — 
400 cm"'^, typical of diffuse molecular clouds, not dense 
clouds. Also, we can estimate the fractional abundance 
of CO (x(CO) = iV(CO)/A^H) in the observed sight 
lines and compare it to the sola r system abundance of 
carbon (x(Ctot) = 2.9 x lO"'*; iLoddersI [200l to de- 
termine if CO is the dominant carbon bearing species. 
We estimate A^h from the color excess (see Table S]), 
and use observed relationshi ps between A^ ( CH), A^(CN), 
and A^(CO) (Sonncntrucker_elalJ [20071 : ISheffer et all 
I2008D. in concert with CH and CN column densities 
(IHirschauer et al.|[200ft) to estimate A^(CO). In the ALS 
8828 and HD 254577 sight lines a;(CO) - 5 x 10"^ and 
2 X 10~6, respectively, much smaller than the assumed 
total carbon budget. In the other 3 sight lines x(CO) 
ranges from about 1 x 10~^ to about 1 x 10"^. These 
estimates show that most carbon is not in the form of 
CO, but does not rule out C as the dominant carbon 
bearing species. To do so, we use the observed relation- 
ship between CO/H2 an d (C-f CO) /Ctot shown in Figure 
6 of lBurgh et all pMl . For CO/H2 - 8 x 10-^— the 
largest value estimated along any of our sight lines — 
observations show that both CO and C account for only 
a small fraction of the total carbon budget, thus indicat- 
ing that carbon is predominantly in ionized form. 

To improve upon these rough arguments though, obser- 
vations yielding the relative abundances of C^, C, and 
CO are necessary. The v—1-0 fundamental and v=2- 
overtone rovibrational bands of CO near 4.6 /im and 
2.3 /im, respectively, can be observed with NIRSPEC 
and IRCS. Various electronic transitions of CO and C i 
are available in the far ultraviolet (1100 A-1700 A), and 
can be observed with either COS or STIS aboard Hubble. 
Finally, a weak intersystem line of C 11 is at 2325 A, and 
may also be observable with COS and/or STIS. Com- 
bined, these observations would allow us to determine 
the predominant carbon-bearing species along each sight 
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line, and give us a better understanding of cloud condi- 
tions being probed. 

6.2. Higher Ionization Rate 

If the gas conditions in all of our observed sight lines 
are similar, then the cosmic-ray ionization rate must be 
higher toward ALS 8828 and HD 254577. Such varied 
ionization rates can be the result of differing cosmic-ray 
fluxes in each sight line. If we assume that the SNR accel- 
erates particles isotropically (i.e. the spectrum of cosmic 
rays leaving the remnant is identical everywhere along 
the blast wave) then the different cosmic-ray spectra op- 
erating in each sight line must be due to propagation 
effects. 

To determine whether or not cosmic rays accelerated 
by IC 443 can even produce the high i nferred ioniza- 
tion ra tes, we use the methods described in llndriolo et al.l 
(|2009[) to compute the expected ionization rate fo r vari- 
ous co smic-r ay spectra. .A bdo ct al. (2010), Torrc s~t al.l 
(HODS'), and iTorres et all (j20ia ) constrain the proton 
spectrum above ~ 100 MeV near IC 443 from the 
observed gamma-ray spectru m. Although the broken 
power-law proton spectrum in lAbdo et al.l (|2010[ ) is given 
as a power law in kinetic energy (flux oc E^^^^^ when 
Ekin < 69 GeV), we change this to a power law in mo- 
mentum (flux oc where pc = [(E'fcm + rUpC^Y — 
{nipC^Yf'^) to account for the fact that diffusive shock 
acceleration is expected to produce a spectrum of this 
form. This s ubstitution only d iffers from the relation 
considered by lAbdo et al.l (|2010t ) in the non-relativistic 
regime where p oc E'^-^. As we extrapolate the spec- 
trum to lower energies then, the particle distribution is 
flattened relative to a pure power law in kinetic energy. 

Integrating this extrapolated spectrum to a low-energy 
cutoff of 5 MeV, we find C2 ~ lO^^'^ s^^. A 5 MeV cutoff 
was used because partic les of this energy hav e a range of 
a few times 10^^ cm~^ (jPadovani et al.|[2009D . similar to 
the sight lines we consider here. Using the d = 10 pc con- 
tinuous injection spe ctrum (also extrapo lated to low en- 
ergies as above) from lTorres et all (j20 08) produces simi- 
lar results, while their ci = 30 pc spectrum fails to repro- 
duce even the ionization rate predicted by the lo cal inter- 
stella r cosmic-ray spectrum ((^2 ~ 4 x 10"^'' s~^: lWebbeH 
|1998() . Note that the difference between the 10 pc and 
30 pc theoretical spectra is noi the result of energy losses, 
but due to the fact that lower energy particles have not 
yet had sufficient time to travel far from IC 443 given its 
age of 30,000 yr. Although these spectra are not particu- 
larly well suited for estimating the cosmic-ray ionization 
rate — they are based on observations which depend on 
processes requiring Ekin > 280 MeV and so are not well 
constrained at energies of a few MeV where ionization 
is much more efficient — they do suggest that cosmic-rays 
accelerated by IC 443 are capable of generating the ion- 
ization rate inferred from Hj^, and also provide indepen- 
dent constraints on the flux of high-energy cosmic rays 
near IC 443 which complement the low-energy compo- 
nent studied in this paper. 

The propagation effect s included in t he model co smic- 
ray spectra presented in iTorres et al.l ()2008l I2010D may 
also be able to explain the differences inferred in (2- 
Those authors suggest that the difference in centroid po- 
sition between the GeV {EGRET and Fermi LAT) and 



TeV {MAGIC and VERITAS) gamma-ray sources can be 
explained by having the gamma rays of different energies 
originate in separate clouds. They propose that the lower 
energy gamma rays arise from tt" decay in a cloud 3-6 pc 
away from the expanding shell of IC 443, and the higher 
energy gamma rays in a cloud about 10 pc in front of the 
SNR. This explanation requires that cosmic-ray propaga- 
tion is energy dependent, such that high energy particles 
have diffused farther away from the SNR than low en- 
ergy particles. In such a model, the cosmic-ray spectrum 
varies as a function of position, and so the ionization rate 
must as well. 

Because lower energy cosmic rays have yet to prop- 
agate very far from IC 443, the ionization rate should 
decrease with increased distance away from the SNR. 
If the clouds probed by ALS 8828 and HD 254577 are 
closer to IC 443 than the clouds probed by the other 
3 sight lines, then the difference in inferred ionization 
rates is easily explained. The positions of our target 
sight lines with respect to IC 443 are shown in Figure 
[11 and the on-sky distances from the center of subshell 
A to each of the background stars is listed in Table 
m Of the 5 sight lines, only HD 43703 is a consider- 
able distance away from the remnant, so differences in 
the remaining 4 sight lines must be due to line-of-sight 
distances. Gas ve locities for the dominan t CH compo- 
nents reported bv iHirschauer et all (|2009D vary by only 
about 3 km s~^ between all of our sight lines, sug- 
gesting that the absorption may arise from the same 
cloud complex, but because IC 443 is located near the 
Galactic anti-center such an analysis is highly uncer- 
tain. The HD 254577 sight line passes through regions of 
HCO+ emissi on ([^ic kman et al. 199j) and H2 em ission 
(jBurton et al.lll988riInoue et al.lll993t iRho et al.l lMT). 
both of which trace shocked gas, and is in close proxim- 
ity to an OH (1720 MHz) maser which requires sho cked 
gas and a high ionization rate (jHewitt et al.l l2006il) . so 
it is plausible that the observed K'^ absorption arises in 
material very close to the SNR shock. The ALS 8828 
sight line, however, is not coincident with shock tracers, 
so it is unclear at this location how close the foreground 
cloud is to the SNR. Still, given the drastic difference in 
the 10 pc and 30 pc cosmic-ray spectra from lTorres et al.l 
(|200l ( see their Figure 1), the gas probed by the 3 sight 
lines without observed Hj^ would not have to be that 
much farther away than the gas probed by ALS 8828 to 
explain the inferred ionization rates; something on the 
order of 10 pc farther away would suffice. 

Aside from the distance between the site of particle ac- 
celeration and the clouds in question, various other prop- 
agation and acceleration effects could account for the dif- 
ference in inferred ionization rates. Cosmic rays diffuse 
through space as they scatter off of Alfven waves which 
are presumed to be generated by the particles themselves. 
For clouds with higher densities of neutral gas, the damp- 
ing of these waves (via ion-neutral collisions) becomes 
more efficient and the streaming velo city of cosmic rays 
thus increases (jPadoan fc Scaloll2005l ). Instead of diffus- 
ing then, particles will free-stream and spend much less 
time in the cloud (i.e. have fewer chances to ionize am- 
bient material). As a result, regions of low gas density 
should be expected to have higher ionization rates than 
regions of high gas density. Another possibility is that 
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the net flux of cosmic rays into a cloud (due to ioniza- 
tion losses, nuclear interactions, etc. within the cloud) 
sets up an anisotropy that causes the growth of Alfven 
waves in the plasma surrounding the cloud. Lower en- 
ergy particles {Ekin less than a few hundred MeV) scat- 
ter off of these waves and are irn peded from entering 
the cloud (jSkilling fc Strong 119760 . Because the parti- 
cles most efficient at ionizing hydrogen are excluded from 
denser clouds, this effect also predicts a higher ionization 
rate in regions of lower density. However, due to the sim- 
ilar densities reported in Hirschauer et al. (2009) for our 
target sight lines, these effects seem unlikely candidates 
for causing the difference in inferred ionization rates. 

The final effect we consider in attempting to explain 
these variations in the ionization rate is the escape of 
cosmic rays upstream from the shock where diffusive 
shock acceleration occurs (i.e. away from the SNR). 
This sub ject has been the focus of several recent stud- 



ies (e.s. ICaprioh et al.l [2001 I20T0I: iReville et all 12001 



lOhira et al.ll2010D which find that particles can escape in 
the upstream direction, although these tend to be only 
the particles with the highest energies. For the discussion 
above, we have assumed that low-energy cosmic rays have 
escaped from the shock and are diffusing away from the 
SNR. However, if low-energy cosmic rays do not escape, 
but are instead preferentially advected downstream (i.e. 
into the SNR) , then the ionization rate in the post-shock 
gas inside the SNR should be higher than in the gas exte- 
rior to the remnant. The differing ionization rates could 
then be explained if the sight lines toward ALS 8828 and 
HD 254577 probed gas interior to IC 443. As mentioned 
above, the HD 254577 sight line is coincident with var- 
ious shock tracers, as well as an OH (1720 MHz) maser 
(which arises from the post-shock gas inside the SNR). 
Additionally, the velocity of this maser, —6.85 km s^^, 
is consistent with the H3 velocities reported in Table [3l 
making it highly plausible that the H|^ absorption to- 
ward HD 254577 arises from shocked gas inside of IC 
443. Consequently, the inability of low-energy cosmic 
rays to escape from IC 443 provides an alternative to the 
diffusion of particles and differing distances between the 
remnant and gas probed by our sight lines in explaining 
the inferred ionization rates. 



6.3. Implications 

Given either of the cases discussed above (low electron 
density or high ionization rate) , we can comment on the 
flux of low-energy cosmic rays accelerated by SNRs. In 
the case that the exceptional H;^ column densities ob- 
served are due to a lower destruction rate (i.e. lower 
electron density), then the ionization rate near IC 443 
is no higher than already found toward various diffuse 
molecular cloud sight lines. This would indicate that 
the flux of low energy cosmic rays near SNRs is not sub- 
stantially different than in the Galactic ISM, and suggest 
that either SNRs are not the primary accelerators of such 
particles, or that low energy particles have yet to escape 
from IC 443. 

In the case that the inferred ionization rates of a few 
times 10~^^ s~^ are correct, IC 443 must be accelerating 
a large population of low-energy cosmic rays. Either this 
population must be escaping upstream from the site of 
diffusive shock acceleration (i.e. traveling outward from 



the SNR shock) such that the clouds closest to the rem- 
nant are experiencing a large flux of cosmic rays, or the 
2 sight lines with H3 detections probe gas inside of IC 
443 where low-energy cosmic rays have been advected 
downstream. In either situation, it is unclear if such a 
population of cosmic rays accelerated by all SNRs within 
the Galaxy will propagate far enough from their sources 
to affect the flux of cosmic rays at some arbitrary posi- 
tion. As a result, it is difhcult to definitively say whether 
or not SNRs are responsible for accelerating the large 
flux of low-energy cosmic rays necessary to produce the 
C2 4 X 10~^^ s~ ionization rate inferred in many diffuse 
Galactic sight lines. 

7. CONCLUSIONS 

We have searched for absorption along 6 sight lines 
that pass through molecular material in the vicinity of 
the SNR IC 443. Two of the observed sight lines, ALS 
8828 and HD 254577, have large column densities of 
, while the other 4 show no absorption features. The 
cosmic-ray ionization rates inferred from the 2 detections 
are a few times lO"""^^ s~^, higher than ever previously 
found in diffuse molecular clouds. Upper limits to the 
ionization rate in the other sight lines, however, are con- 
sistent with values found along average Galactic sight 
lines of about 4 x 10^^^ s^^. These differences may be 
due to overestimates of the electron fraction or cosmic- 
ray propagation and acceleration effects, but the com- 
plexity of the region makes it difficult to attribute the 
results to any one cause. Future observations of C 11, 
C I, and CO toward our target sight lines should allow 
us to better discriminate between the two possibilities, 
and thus determine whether or not IC 443 produces a 
large flux of low-energy cosmic rays. 

In addition, surveys of H;^ near IC 443 and other SNRs 
thought to be interacting with molecular clouds (e.g. 
Vela, W 28, W 44, W 51C) should allow us to further in- 
vestigate cosmic-ray acceleration in such environments. 
By more extensively mapping H;^ absorption near super- 
nova remnants, we can determine where the H^ resides 
(interior post-shock gas or exterior pre-shock gas), and 
thus where the flux of low-energy cosmic rays is high- 
est. Such observations may also provide insight into the 
efficiency with which accelerated particles are advected 
downstream into remnants, and so add important con- 
straints to models of cosmic-ray acceleration. 
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Fig. 1. — This image of IC 443 is from the Second Palomar Observatory Sky Survey (POSS-II) using the red filter, and was obtained from 
the STScI Digitized Sky Survey. Target background stars are to the immediate left of the uppercase letters, and are labeled as follows; 
A-ALS 8828; B-HD 254577; C-HD 254755; D-HD 43582; E-HD 43703; F-HD 43907. The cross marks the center of subshell A (the shell 
on the northeast side of IC 443) at a = OG^lT^OS.f , S = +22°36'39.4" J2000.0. The three black squares mark the positions of OH 
(1720 MHz) maser emission reported bv lHewitt et al.l l|200g ). 
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Fig. 2. — Model atmospheric transmission spectra at wavelengths near targeted Hg" transitions, generated by ATRAN I lLordI [T992I ) 
assuming observations were performed from Mauna Kea, 1.6 mm precipitable water vapor, and an air mass of 1.15, and smoothed to a 
resolving power of about 25,000. Dashed vertical lines mark the rest positions of Hg" transitions, and are labeled accordingly. The bottom 
axes give the rest wavelengths in ^m, while the top axes show the effect of relative motion between the observer and interstellar gas. Average 
velocity shifts (due to the Earth's motion and interstellar gas motion) for the spectra taken at Keck and Subaru are about —18 km s~^ 
and —2 km s~^, respectively. 
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Fig. 3. — Spectra of six stars in the IC 443 region covering various Hg transitions. Vertical dashed lines show the expected position of 
absorption lines due to Hg' given the velocities reported by Hirsch auer et all l|2009) for cloud components with the most CH absorption. 
The shorter dashed line in the R{1, 1)"--R{1, 0) spectra shows the position of the R{1, 1)" line, which is 36 km s"'^ away from the R{1, 0) 
line. For HD 254755, HD 43582, HD 43703, and HD 254577, the _R{1, 1)"--R(1, 0) spectra are the combination of NIRSPEC and IRCS data, 
while the _R(1, 1)', Q(l, 1), and Q(l, 0) spectra are only from IRCS. All of the ALS 8828 spectra are from NIRSPEC, and HD 43907 spectra 
are from IRCS. Of these six sight lines, only HD 254577 and ALS 8828 show H;^ absorption features. The _R(1, 1)", R{1, 0), 1)' , and 
Q(1,0) lines are visible toward HD 254577 (the lower S/N in the Q{1, 1) spectrum and smaller dipole moment of that transition precludes 
its detection). For ALS 8828, relatively strong R{1, 1)" and -R(l, 0) lines are visible. Even with a large amount of Hg" along the sight line 
though, the R{2, 1)" and -R(3, 3)' transitions arising from higher energy states are not detected. 



11 



TABLE 1 

Hg" Transition Properties 



Wavelength \fi\^ 



Transition 


(^tm) 


(D2) 




3.668083 


0.0158 


R{1,0) 


3.668516 


0.0259 




3.715479 


0.0141 




3.928625 


0.0128 


Q(1,0) 


3.953000 


0.0254 


R{3,3)' 


3.533666 


0.0191 


Note. — 


Wavelengths and dipole 



moments for targeted transitions in 
the 1/2^0 band of Hg" (values 
from lGoto et al.l2002l . and references 
therein). 
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TABLE 2 
Observations 



Object 


Date(s) of Observation 


Telescope 


Integration Time 
(min) 


ALb 8828 


2009 Nov 5 


Keck 


40 




2009 JNov D 


Keck 


24 


riU 254577 


2009 JNov 5 


Keck 


40 




2009 Nov 6 


Keck 


20 




2009 Dec 12 


Subaru 


96 




2009 Dec 13 


Subaru 


36 


HD 254755 


2009 Nov 5 


Keck 


40 




2009 Nov 6 


Keck 


20 




2009 Dec 13 


Subaru 


120 


HD 43582 


2009 Nov 5 


Keck 


28 




2009 Nov 6 


Keck 


40 




2009 Dec 12 


Subaru 


120 


HD 43703 


2009 Nov 6 


Keck 


68 




2009 Dec 13 


Subaru 


120 


HD 43907 


2009 Dec 12 


Subaru 


42 
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TABLE 3 
Absorption Line Parameters 







''LSR 


b WHM 


Wx 


a(Wx) 


N(J, K) 


a{N) 


Object 


Transition 


(km s ) 


(km s ) 


(10~^ tim) 


(10^*' fim) 


(10 cm ^) 


(10 cm ^) 


ALS 8828 


_R(1, 1)" 


-6.0 


14.0 


6.4 


0.5 


2.64 


0.22 




-R(l, 0) 


-8.3 


14.1 


7.0 


0.5 


1.76 


0.14 


HD 254577 


_R(1, 1)" 


-6.7 


24.9 


3.0 


0.4 


1.22 


0.15 




R{1, 0) 


-6.2 


16.2 


3.8 


0.3 


0.95 


0.07 




-R(l. iV 


-6.2 


27.0 


3.4 


0.4 


1.55 


0.17 








16 


< 3.2 




< 1.54 






Q(l, 0) 


-6.1 


14.3 


2.5 


0.5 


0.59 


0.12 


HD 254755 


_R(1, 1)" 




16 


< 1.2 




< 0.49 






R(l, 0) 




16 


< 1.2 




< 0.30 






-R(l, iV 




16 


< 1.4 




< 0.64 






OCl 1 ~) 




16 


< 3.3 




< 1.58 






oh 0) 




16 


< 1.5 




< 0.36 




HD 43582 


_R(1, 1)" 




16 


< 1.9 




< 0.78 






R{1,0) 




16 


< 1.9 




< 0.48 






1)' 




16 


< 1.5 




< 0.68 






<9(1,1) 




16 


< 3.2 




< 1.55 






Q(1,0) 




16 


< 3.0 




< 0.71 




HD 43703 


R(l,l)" 




16 


< 1.2 




< 0.50 






R{1,0) 




16 


< 1.2 




< 0.31 






R{1, ly 




16 


< 1.3 




< 0.62 






Q(l,l) 




16 


< 3.5 




< 1.68 






Q(1,0) 




16 


< 2.3 




< 0.56 




HD 43907 


_R(1, 1)" 




16 


< 4.4 




< 1.83 






R{1,0) 




16 


< 4.4 




< 1.12 






R{1, ly 




16 


< 5.3 




< 2.41 





Note. — Coluirm 3 (t'LSR ) gives the interstellar gas velocity in the local standard of rest frame. Colunm 
4 (FWHM) gives the full width at half-maximum of tfie absorption features. In tlic case of non-detections, 
the FWHM was set to 16 km s^^, the resolving power of IRCS on Subaru in our particular setup, for the 
purpose of computing column density upper limits. Columns 5 and 6 show the equivalent width, Wx, and 
its la uncertainty, a(Wx), respectively. Upper limits to Wx are equal to 3cr{Wx)- Columns 7 and 8 give the 
column density of H^ in the state each transition probes, N{J, K), and its uncertainty, a-{N), respectively. 
Upper limits to the Hg" column density are equal to 3a{N). 
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TABLE 4 
Target Sight Line Properties 



Target 


r 
(pc) 


T 
(K) 


(cm-'^) 


(lO^l cm- 


L 

'") (pc) 


Af{H2) 
(lO^l cm-2) 


^(Hj) 
(1014 cm-2) 


C2 

(10-'«^ s-') 


ALS 8828 


6.8=^ 


60'' 


300= 


3.0 


3.2 


2 1+'* 


4.4 ±0.26 


16+* 


HD 254577 


7.0^ 


35<i 


325'' 


3.6 


3.6 




2.2 ±0.34 




HD 254755 


8.6 


35'^ 


200'' 


2.5 


4.0 


1 1+°-^ 


< 0.6 


< 3.5 


HD 43582 


5.4=^ 


60*' 


200" 


1.9 


3.0 


o.5«:| 


< 0.8 


< 9.0 


HD 43703 


13.9 


60'' 


300^^ 


1.7 


1.8 


0.8l°J 


< 0.6 


< 5.7 


HD 43907 


20.6 


60'' 


300^^ 


1.5 


1.6 


o.4i°:3 


< 2.1 


< 40 


Note. — Various parameti 


ars used 


in our analysis for the target sight hues 


in this study. Column 2 gives 



the on-sky distance (r) from the center of shell A of IC 443 (a = 06'^17™08.4=, 5 = +22°36'39.4" J2000.0) 
to each sight line assuming the remnant is at a distance of 1.5 kpc. The radius of shell A is about 7 pc. 
Temperatures (T) were taken from the C2 rotation-excitation analysis when available or set to 60 K, and 
number densities (nn) were taken from either the C2 rotation-excitation analysis or restricted chemical 
analysis, both reported in Hirschauer et aP ( 2009). Uncertainties in Jipj are taken to be ±100 cm""^. Color 
excesses were decreased by 0.3 mag to remove the contribution from foreground gas (color excesses and the 
foreground correction are given in Hirschauer et al. 20Q3), and the relationship A'^h ~ i?(-B— V) x5.8x 10^' 
cm" mag-' IjBohlin et al.lll978l : iRachford et al.li2003 l was used to compute the total hydrogen column 
densities, A^h- Path lengths (L) were calculated from njj and Ajj. Molecular hydrogen column densities 
were calculated from the relationship A(CH)/A(H2) = 3.5ti;4 X IQ-* l lSheffer et al.ll200g) using the 
dominant CH components from Hir schauer et al] I I2009I1 . Uncertainties in Af(H2) are dominated by the 
scatter in the above relationship, not uncertainties in A(CH). Upper limits for Af(Hg') are the 3(t 
uncertainties from the observations, and the upper limits for ^2 are based solely on those values. To 
account for the uncertainty in Af(H2), the upper limits for 1^2 should be multiplied by 1.5. 

^ This sight lines passes through the remnant. 

'' Kinetic temperature assumed in restricted chemical analysis. 

Average density from restricted chemical analysis. 
'' Determined from C2 rotation-excitation analysis. 

As no CN was detected toward HD 43582, here we have taken the lower bound on the density from the 
restricted chemical analysis for other sight lines. 



